Magnetic response of the J1 — J2 spin Hamiltonian from 
classical Monte Carlo and Schwinger boson mean field theory 
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Magnetic susceptibilities at several potential ordering wave vectors (0,0), (vr,0), and (tTjTt) are 
analyzed for the antiferromagnetic J1—J2 spin Hamiltonian by classical Monte Carlo and Schwinger 
boson mean field theories over the parameter range < 2J2/J1 < 2. We find a nearly linear-T 
behavior of the uniform susceptibility that extends up to the temperature scale of Ji within both 
calculation schemes when 2J2/J1 is sufficiently removed from the critical point 2J2/J1 = f. The 
window of linear temperature dependence diminishes as the critical point is approached. 
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I. INTRODUCTION 

One of the outstanding issues for the spin density 
wave (SDW) phase of the parent FeAs compound is 
whether a localized spin model can be justifiably used 
for the description of spin dynamics. The concern arises 
mainly from the fact that, unlike the parent compound 
of cuprate high-temperature superconductors, the par- 
ent material for FeAs-based superconductors is metallic 
and the validity of the description which ignores the itin- 
erant electrons is not nearly as water-tight as with the 
half-filled cuprate compound. Despite these concerns, 
the Ji — J2 spin model had been proposed for the spin 
dynamics of the FeAs parent compound^'^ which, dis- 
regarding the multi-orbital character of the FeAs band 
structure, reads 



Hs — Ji ^ Si 

iv) 



(1) 



The nearest neighbor (NN) and next-nearest neighbor 
(NNN) pairs are written as (ij) and (ik), respectively. 
From first-principles calculations one obtains the esti- 
mate Ji ~ J2 and Ji in the range of several hundred 
Kelvin^^'^. 

Of particular relevance in connection with the mag- 
netic response of FeAs is the behavior of the uniform 
susceptibility Xu- The temperature dependence of Xu 
was found to be Hnear up to 500-700K'*, well beyond 
the SDW ordering temperature Tsdw ~ 150K. The rel- 
evance of the Ji — J2 spin Hamiltonian to the FeAs 
phenomenology largely depends on its ability to re- 
produce the observed temperature dependence of x«- 
The uniform susceptibility for the Ji-only model has 
been calculated semi-classically by Takahashi^ and more 
recently by classical Monte Carlo (MC) calculations^, 
demonstrating a roughly linear behavior of Xu{T) in the 
low-tcmperature regime below the mean-field transition 
temperature, Tmf- The quantum spin-1/2 Heisenberg 
Hamiltonian was also found to have a low-temperature 



region with a roughly linear temperature-dependent Xu' ■ 

A recent note by Zhang et al.^ cites such susceptibil- 
ity behavior as evidence of pre-formed moments in FeAs 
in the paramagnetic regime. A quasi-linear tempera- 
ture dependence of Xu in the Ji — J2 spin Hamiltonian 
was demonstrated in a spin-wave calculation based on 
the Dyson-Maleev transformation of the spin operators*. 
A similar temperature dependence of the susceptibility 
was also found in a model based on the nonlinear sigma 
model description of the spin dynamics^. Both ideas*'^ 
support, and are based on the presence of pre-formed 
magnetic moments in the paramagnetic phase^'^. In con- 
trast to the Curie- Weiss susceptibility behavior of a local 
moment, the linear-T susceptibility is a reflection of the 
non-trivial correlation among the moments in the para- 
magnetic state, and a deeper understanding of the mag- 
netic response in such a correlated paramagnetic phase 
is highly desirable. 

A lot is known about the classical Ji — J2 model from 
previous theoretical works^°'^^. Although no true long- 
range order is possible in two dimensions, the mean-field 
ground state is antiferromagnetically ordered at the wave 
vector k = (7r,7r) when 23^1 J\ < 1. For 2J2/J1 > 1, 
the two sublattices separately support a mean-field stag- 
gered magnetic phase while the relative orientation of 
the sublattice magnetization vectors remains undeter- 
mined. It is the thermal fluctuations which selects a 
coUinear spin order at fc = (tt, 0) or fc = (0, tt) through 
a order-by-disorder mechanism^^. A quantum spin-1/2 
version of the Ji — J2 has recently received some the- 
oretical attention^^, but the relevant uniform magnetic 
susceptibility Xu has not been calculated in the quantum 
model. 

In this paper, we re-visit the susceptibility behavior of 
the J1—J2 spin model in a comprehensive and unified fash- 
ion using both classical Monte Carlo (MC) methods and 
the Schwinger boson mean field theory (SBMFT). Previ- 
ous analytic and numerical works^^*'^^ typically have ex- 
amined the susceptibility for a particular ordering wave 
vector. In the J1 — J2 spin model the dominant ordering 
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wave vector switches from (tt, tt) for 2J2/J1 < 1 to (7r,0) 
for 2J2/J1 > 1. The fate of the non-dominant suscep- 
tibihtics, for instance the susceptibihty at fc = (tt, 0) for 
the 2J2/J1 < 1 region, has not been addressed carefully. 
The uniform susceptibility x« is not expected to be dom- 
inant for any 2J2/J1 ratio, but its behavior has not been 
understood in a systematic way, either. 

At the moment the controversy surrounding the 
itinerant^^ vs. localized moment origin®'^ of the linear-T 
susceptibility is an on-going issue. A careful and system- 
atic examination of the magnetic response of the Ji — J2 
spin model by employing well-established techniques is 
expected to be a helpful addition to our knowledge re- 
quired to settle the issue in the long run. The orga- 
nization of the paper is as follows. In Sec. II, Monte 
Carlo calculation of the magnetic susceptibilities at sev- 
eral wave vectors are presented for all IJij J\ ratios be- 
tween and 2, which encompasses the Hcisenberg limit 
2J2/J1 =0 and the critical point at 2J2/ J\ = 1. And in 
Sec. Ill a Schwinger boson mean-field calculation employ- 
ing several order parameters which embody both (tt, tt) 
and (77,0) ordering are presented over the same 2J2/J1 
range. We conclude with a discussion in Sec. IV. The 
properties of the uniform and (tt, 7r)-staggered magnetic 
susceptibilities of the Ji-only Hamiltonian in the SBMFT 
are reviewed in the Appendix. 



II. MONTE CARLO EVALUATION OF THE 
MAGNETIC SUSCEPTIBILITY 

Classical mean field analysis of the J1—J2 model assum- 
ing the ordering pattern {Si) ~ mi(— 1)^' -|-m2(— l)^*"*"^' 
(coexistence of (tt, 0) and (tt, tt) order) is reduced to the 
coupled set of mean-field equations 



mi+m2 = tanh|^— [(J2mi-|-(Ji-J2)m2]J , 
mi-m2 = tanh^^[J2mi-(Ji-J2)m2]^ . (2) 

As shown in Fig. 1, either a (tt, 7r)-ordered phase (mi = 
0,m2 ^ 0) or a (tt, 0)-ordered phase {mi 7^ 0,m2 = 0) 
is favored depending on the 2J2/J1 value being less or 
greater than one. No co-existence region exists where 
mi and m2 are simultaneously nonzero. A first-order line 
separates the two phases at 2 J2/ Ji = 1. In the Ginzburg- 
Landau picture, a repulsive interaction u|mip|m2p term 
with V > exists between the order parameters. 

With this classical mean field phase diagram in mind, 
we have measured the magnetic susceptibilities at three 
wave vector positions, k ~ (0, 0), (tt, 0), and (tt, tt), for 
the range < 2J2/J1 < 2 that also encompasses the 
critical point 2J2/J1 = 1. For a given k vector, the sus- 
ceptibility Xfc for the component a = a;, y, 2; is computed 
according to 




2J2/J1 



FIG. 1: Mean-field phase diagram of the J1 — J2 classical 
spin model, (tt, 7r)-ordered phase is separated by a first-order 

line (dotted) at 2J2/J1 = 1 from the (tt, 0)-ordered phase. 
Crossover temperature Tco extracted from the uniform sus- 
ceptibility in the classical Monte Carlo calculation is shown 
for comparison. 



i 

Here N stands for the number of lattice sites, and (• • • ) 
is the thermal average. 

In agreement with the mean-field analysis, the (tt, tt) 
susceptibility rises sharply at a finite temperature Tc in 
the region 2J2/J1 < 1. For the (7r,0) susceptibility the 
sharp rise occurs on the 2J2/J1 > 1 side. At 2J2/J1 = 1 
both susceptibilities showed increasing behavior at lower 
temperatures. This temperature however decreases 
when a larger lattice size is used in the simulation, in ac- 
cordance with the Mermin- Wagner theorem which pro- 
hibits long-range order in a continuous spin model. We 
will not discuss the behavior of the dominant susceptibil- 
ities any further and turn our attention to non-dominant 
susceptibilities. 

Results of the MC calculation for the non-divergent 
susceptibilities are shown in Fig. 2. The uniform suscep- 
tibility X(o,o)j which remains non-divergent for all val- 
ues of 2J2/J1, does show the low-temperature linear-T 
behavior up to a certain crossover temperature scale, 
Tco^^- Furthermore, Tco grows as the 2J2 — J1 value 
moves away from criticality. The crossover temperature 
Tco c&n be extracted from the X(o.o) data due to the 
fairly abrupt deviation from low-temperature linearity at 
T w Tco- A plot of Tco over the range < 2J2/J1 < 2 
is shown in Fig. 1 along with the mean-field transition 
temperature Tmf- Due to the thermal fluctuation ef- 
fects, Tco remains substantially smaller than the corre- 
sponding Tmf- Even at the critical point 2J2/J1 = 1 the 
crossover scale remains finite at Tqo/Ji '^0.1. 

The behavior of non-divergent X{w,o) for 2J2/J1 < 1, 
and X(7r,7r) for 2J2/J1 > 1 are also shown in Fig. 2. In- 
terestingly, these susceptibilities also show a linear tem- 
perature dependence at low temperatures, and find the 
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(a) 2J2/J1 < 1 



(b) 2J2/J1 > 1 




(e) Slope of Zu 



(f) Slope of X(i,,o) and xf,,,) 
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FIG. 2: (color online) Temperature dependence of (a) uniform 
susceptibility xfo.o) 2J2/J1 < 1, (b) X{o,o) 2J2/J1 > 1, 
(c) staggered susceptibility X^tt.o) ^^o'' 2J2/J1 < 1, and (d) 
staggered susceptibility X(tt -k) f^ir 2J2/J1 > 1. Data are ob- 
tained for a 16x16 lattice with 5 x 10^ MC steps per data 
point. Several calculations run on a 40x40 lattice gave vir- 
tually identical curves as those of a 16x16 lattice over the 
whole temperature interval. The crossover temperature such 
as indicated by an arrow in (a) and (b) marks the devia- 
tion from low-temperature linearity. We have plotted the z- 
component of the susceptibilities for an easier determination 
of Tco^'^- (e) Slopes of the low-temperature linear part of 
X(o 0) vs. T/Ji. (f) Slopes of the low-temperature linear part 
of X(^,o) for 2J2/J1 < 1 and X{n,-^) for 2J2/J1 > 1. Both 
slopes become sharper near the critical point. 



maximum value at T rs Tqo ■ Unlike the uniform suscep- 
tibility, these two quantities exhibit divergent behavior 
right at the critical point 2J2/J1 — 1, and are omitted 
from the plot in Fig. 2. In terms of absolute scale, the 
(tt, 0) and (tt, tt) susceptibilities are about an order of 
magnitude bigger than Xu- The linearity can be char- 
acterized by fitting the low-temperature part with the 
form = B(T/Ji). The slope B of the linear parts of 
the susceptibilities are shown in Fig. 2(e) for the uniform 
part and Fig. 2 for the staggered parts. 



Based on the results of MC calculations, we conclude 
that the linear-T behavior of the uniform susceptibility is 
a natural aspect of the Ji — J2 spin model for the 2J2/J1 
ratio that spans both sides of the critical point, and pre- 
sumably for a wider class of spin models whose mean- field 
ordering wave vector occurs well away from (0,0). In 
fact, the linearity of XuiT) persists over a wider temper- 
ature range when 2 J2/ Ji departs further from the critical 
value. If we take the estimate Ji = J2 = 500K^, the ex- 
pected crossover temperature from the MC calculation is 
Too ~ O.8J1 - 400K, in fair proximity to the 500-700 K 
up to which the linear susceptibility has been observed. 
It is also clear that the linear-T susceptibility behavior 
alone does not imply the proximity of the FeAs system 
to a critical point nor the existence of high degree of spin 
frustration. 



III. SCHWINGER BOSON MEAN FIELD 
THEORY OF THE MAGNETIC 
SUSCEPTIBILITY 



The quantum analogue of the classical susceptibilities 
can be obtained by treating the Ji — J2 spin Hamilto- 
nian in the framework of Schwinger boson mean field the- 
ory (SBMFT). A prior calculation* assumed an ordering 
wave vector = (tt, 0), or equivalently at Ky = (0,7r), 
and ignored the possibility of ordering at Q — (7r,7r). 
While this assumption may be sound deep in the J2- 
dominated part of the phase diagram, here we want to 
make an unbiased statement that spans both small and 
large 2J2/J1 regions and keep both (tt, 0) and (7r,7r) or- 
der parameters in the theory. The details of the self- 
consistent Schwinger boson calculation for the Ji-only 
spin model can be found in the work of Auerbach and 
Arovas^'"^. 
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FIG. 3: (color online) Schematic picture of the pairing ampli- 
tudes Aij used in our calculation. The values of Ai-, for the 
arrows going from site i to site j are indicated on the right. 
Note that A^i = — Ai, . 



For SBMFT calculations one first re- writes the J\ - 
Hamiltonian in the mean-field form 
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MF 



+ 



(4) 



where the operators are Ay = 6ii6j2 — ^i2&ji = —-^ji, 
and Ajj = {Aij). The 6,i and 6,2 are the two species 
of boson operators that obey the occupation constraint 
feji^ii + 6^26,2 = 25', which is enforced with the Lagrange 
multipUer Aj. 

The distribution of Ay adopted in our calculation is 
schematically shown in Fig. 3. The nearest-neighbor 
pairing fields are assumed to be modulated at both 
and Q with the respective amplitudes Ai and A2: 



Am±£ = A2(-ir+^' + Ai,(-ir% 

A,,,±^ = A2(-l)^-'+^"^ + Al^(-l)-^ (5) 

By symmetry, it turns out, Aiy vanishes identically and 
we can re- write Ai^^ = Ai. For the assumed ordering at 
(0, tt) wc would have Ai^ = instead. Without loss of 
generality we choose to work with the (tt, 0) modulation 
below. As for the interactions among the next-nearest 
neighbors let us write 



-A, 



(-ir- (6) 



It will be shown that A^ = As whenever they are 
nonzero. 

Due to the enlarged unit cell in real space assumed 
by the ordering pattern, the Brillouin zone is reduced to 
[— ■|,-|]x[— ^,|], and the Hamiltonian is expressed with 
an 8-component spinor 



/ bki \ 

bk+K^ 1 

bk+Ky 1 

bk-{Qi 



fe2 
bl- 

bi 

k+Ky2 



H 



1. \ 



A 



i>k. (7) 



The 4x4 matrix Tik is given by 



/ 

Aife+Aa/; 

, ^Zk 
\ A2k 



where 



Aifc — Asfe — Ag^ 


-A2fc 


A^, 


A'zk 


-Afk 


-Au+A 


"Ag^, Alfe-Asfe 






3fe 



,(8) 



Aife 

A2fe 

A2k 

Aafe 
A' 



3k 



JiAi cos kx, 

Ji A2 (cos kx + COS ky ) , 
Ji A2 (cos kx — cos ky ) , 
J2{Aa — Ab) cos kx cos ky, 
J2(A^+Ab) cos kx cos ky. 



(9) 

Note the antisymmetry {'Hk)fiL> = —{'Hk)vn- The sum 
'}2k extends over the reduced Brillouin zone (RBZ) only. 
The four kinds of momentum "flavors" (fc, + ii'a;, A: + 
Ky, k + Q) can be organized with an index a = 1 ~ 

4: {bkcr,bk+K^a, bk+Kya, bk-iQ(T)^{bk,la,bk,2a, bk,3a,bk,4a)- 

Rotation to the eigenoperators are implemented through 



bk,al = ^(uk.aa-YkM + '"l,aan'k,a2) ' 

a=l 
4 

h,a2 = 5I("'=.««'^M2 -<aa7fe;al), (10) 

which would yield 



H = J2 ^k,a (lk,al'yk,al+^ta2^k,a2) ■ (H) 

kja 



It is convenient to group the wave function compo- 
nents as Ufc,a = (Mfe,la,Ufc,2a,Ufc,3a,Ufc,4a)'^, and Vk,a = 

{vkSa,Vk,2a,Vk^3a,Vk,4,aY' , after which the equations of 
motion becomes 



AUfc,o + "^fe Vfe^a = Ek,aUk,a, 

AVfe,o -I- HkUk,a = -Ek,a^k,a, (12) 



or 



(A^— Sfc o)Ufe^a = 'H^'HkUk,a, 
{X'-EljVk,a = HkH+Vk^a. (13) 

The matrix A4k = TC^'Hk = 'Hk'H'^ is Hcrmitian and all 
its eigenvalues arc non-negative. It proved to be an enor- 
mous advantage to deal with the ordinary Hermitian ma- 
trix problem of Eq. (13) over a non-Hermitian one posed 
by Eq. (12). With the eigenvalues of Mk written m| ^j, 

the energy spectrum is obtained as Ek,a = \jx'^ — w| ^. 
One can either solve for « or Vk a 

from Eq. (13), 

and derive the other part using 



Vfe.a = - 



'HkUk^a 



or Uk,a = 



The normalization for Uk,a and v^^a reads 



(14) 
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"fc.a-"M = nJT ^"i" Vfc -Vfc^b = ' ^ab, (15) 



2Ek,a 



to ensure ,j • Ufe,6-v^^„ • Vfe,fe = 5a6- 

The averages of the boson operators we need are 



where Bk^a is the Bosc-Einstcin function of energy Ek.a, 
and k = —k. The order parameters introduced for the 
mean field calculation are 



4 

(l + Bk.a 



a=l 



(16) 



25+1 



PEk^a 



2 ^ 

k 

= ;^Xl^°^'^^('"'='l«'^^.2a~'"'=.2a^Ma + "fe,3at^fe,4a-'"fc,4a'?^fc,3a)c0th ( 

k,a ^ 

A2 = -^Y^{cosk^+cosky){hkihT^_^2) 

k 

I ' 

^ AtXI [(cOsA;a; + COsA;y)(Ufe,iaVfe_4„-Ufe,4at^fc,la) + (cOsfcx-COsA;y)(Ufc,3a^^fe,2a-'"fc,2at^fe,3a)] COth 



k,a 



l3Ek,a 



= -Jir }_^COSk^COSky{bklb^+K^2 + '^klh+K^2) 



N 

k 

2 . ^ / l3Ek 

— Y^COSk^ cos ky[{Uk,la + Uk,4a)(.^k,2a + V*k,3a) - iUk,2a + Uk,3a){v*k,la + '"k,4a)] COt^ 



Ab = 



N 



— ^ cos fcx COS ky [{Uk,2a - Uk,ia){v*k^la " ^^Ma) " ("fe.la " Wfe,4a) (t^fe,2a -^^fc,3a)] COth 



I ^ X 
= — > -= — coth 



^ ( P^k,a 



N ^ Ek,a V 2 



PEk.a 



r 



(17) 



The last line reflects the boson occupation numbers sat- 
isfy (l/AT) ^,(6t_^6M + &I,2&m) = 25- 

We have considered two spin values, S — 1/2, and 
5 = 1, in the self-consistent calculation. The low- 
temperature (T/Ji = 0.001) behavior of the gap param- 
eters arc shown in Fig. 4 over the range < 2J2/J1 < 2. 
There are two regions showing distinct behaviors, sep- 
arated at 2J2/J1 = 7?c, rjc is a spin-dependent critical 
point. For S = 1/2 the critical point is r]c = 1-2, but 
r]c = 1-1 when S = 1^^. We found that 77c = 1 as long 
as S > 2, which is in agreement with the classical re- 
sult. In Figs. 4 (a) and (b), for 2J2/J1 < rj^ one finds 
/S.A = = As, and Ai 7^ 7^ A2, and 2J2/ Ji > rjc gives 



A^ = Ab nonzero, but Ai = 0. On the other hand A2 
remains nonzero throughout the whole range, although 
its magnitude decreases with the onset of nonzero Aa- 

A surprising feature of the SBMFT result that needs 
to be emphasized is that Ai, representing (77,6) order 
for the NN bonds in the x-direction, becomes nonzero 
even at 2J2/J1 = and stays finite (constant, in fact) all 
through 2J2/J1 = rjc- Their temperature dependence for 
2J2/J1 < r]c is shown in Fig. 4 (c) and (d), which tell us 
that both gap parameters become non-zero at the same 
temperature. A spontaneous loss of rotational symmetry 
is implied by nonzero Ai. As long as A^ = = As the 
Ji — J2 Hamiltonian in SBMFT is the same as the Ji- 
only model, with no infiuence from the presence of the 



6 



1.2 

u 
1 

g 0.6 



5" 0.3 - 

o 





(a)S=l/2 - 


Ai 




Ai 




Aa = Ab 









(b) S=l 


A2 




A, 




Aa = Ab 







0.6 0.9 1.2 1.5 



2.1 0.6 0.8 1.0 1.2 1.4 
2Ji/J, 






(d) S=1,2J2/J,<1.1 


- A2 








• Ai _ 


\ 




i 

1 





0.0 0.3 0.6 0.9 1.2 1.5 0.0 0.3 0.6 0.9 1.2 



T/S(S+1)J, 



T/S(S+1)J, 



FIG. 4: (color online) Dependence of the gap parameters on 
the ratio 2J2/J1 at low temperature, T/Ji — 0.001 for (a) 
5=1/2 and (b) 5 = 1. Temperature dependence of the gap 
parameters Ai and A2 when 2J2/J1 < rjc for (c) 5=1/2 
and (d) 5 = 1. The temperature is divided by the exchange 
energy scale Ji5(5+ 1). The plot shows the onset of both Ai 
and A2 at the same temperature T^^^'^ . 




^ 6 

I" 4 

g. 2 
u 

w 



-2 -1 

/«(T/S(S+l)Ji) 

{c)S=l/2, fa(x„J 




t -3 -2 -1 
/«(T/S(S+l)Ji) 

'(f)S=\, In (x ). 



-3 -2 -1 

/«(T/S(S+1)J,) 




-3 -2-10 1 
;n(T/S(S+l)Ji) 



J2 interaction. This explains the constancy of the gap 
parameters throughout < 2J2/J1 < rjc- 

The equal-time magnetic susceptibility at wave vector 
k given in Eq. (3) can be calculated, where (S^) = and 
thus 



FIG. 5: (color online) Temperature dependence of the suscep- 
tibilities for several 2J2/J1 ratios for 5=1/2 (left column) 
and 5 = 1 (right column), (a) and (d): Xco.o), (b) and (e) 
X{7v,o)- (c) and (f): X(vr,,r)- X{^,o) and X(7v,-^) are shown on the 
log- log scale to emphasize the 1/T divergence at low temper- 
ature. 



(18) 

The SU(2) invariance of the spin-spin correlation was 
assumed in eliminating the {SfS^) correlation from the 

above. The structure factor 5+ = ji^t Sj)e^''<''^^''=^ 
in SBMFT reads 

ki ,k2 

(^LlC2)(^'=i-fc^2fcfe24*,l) + (&Llfefe2,l)(^fel-fe,2&L_fc,2>)- 

(19) 

The first term on the right-hand side, corresponding to 
the spin flip process, is identically zero. The other struc- 
ture factor = J2i ji^r is obtained from 
= S^j,, and the magnetic susceptibility reads 

Xk = ^{s+ + S+,). (20) 



A factor 3/2 has been divided out in arriving at Eq. (20) 
in order to obtain the correct high-temperature results^^. 

Three susceptibilities, X(o,o)j X(ir.o)7 a-nd X(7r,7r)i for spin 
S' = 1/2 and 5=1, were calculated as shown in Fig. 5 
for several 2 J2/ Ji ratios. Above the mean-field transition 
temperature T^^^"^, all the gap parameters vanish iden- 
tically and all the susceptibilities obey the Curie- Weiss 
form Xfe(T) = JiS'(S'-|- l)/r irrespective of k. Due to the 
constancy of the gap parameters for 2J2/J1 < rjc there 
are no variations in the susceptibility in this parameter 
range and it is sufficient to show only one susceptibility 
from this region in Fig. 5. We find T|J|f^^"^ depends 
on the spin S in the manner T^^p^^T _ tJiS{S + 1) 
where i is a 2 J2/ Ji-dependent number. The staggered 
susceptibilities are plotted on a log-log scale in Fig. 5 to 
emphasize the 1/T divergence at low temperatures. 

In the Appendix we show that the uniform suscepti- 
bility in the Ji model follows the form X(o,o) = ^ + BT 
at low temperatures. Clearly the features of the uni- 
form susceptibilities for all ranges of 2J2/J1 are similar, 
with the main difference lying in the values of T^^^'^ 
at which the collapse of the gap parameters occur and 
the susceptibility follows the Curie- Weiss form. The re- 
duction in the slope B (defined in a loose sense as the 
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amount of increase divided by the temperature inter- 
val) with increasing 2J2/J1 is found in both SBMFT 
and MC calculations. Although not nearly linear, the 
monotonic increase of the uniform susceptibility occurs 
over a larger temperature window when 2 J2/ Ji is further 
removed from the critical point, in agreement with the 
earlier findings of Monte Carlo calculation. 

The two staggered susceptibilities X(7r,o) and X(7r,7r) 
show the temperature dependence which are qualita- 
tively similar for all 2 J2/ Ji ranges studied and thus bears 
strong resemblance to that of X(7r,7r) in the Ji model. The 
properties of the latter quantity are discussed in the Ap- 
pendix. Both susceptibilities diverge at low temperature 
as 



X{n,o) ~ ^, and X{^,n) ~ ^ (21) 

on both sides of the critical ratio rjc- The reason for the 
appearance of the factor N is explained in the Appendix 
as well. The reason for the simultaneous divergence of the 
two susceptibilities is the presence of Ai for 2J2/J1 < rjc 
and Aa{b) for 2J2/J1 > rjc, either of which would imply 
(quasi) ordering at k = (tt, 0). Due to the impossibility 
of long range ordering in two dimensions, the associated 
magnetic susceptibilities continue to diverge down to zero 
temperature. The scaling of the susceptibility with A'' is 
tied to the near Bose-Einstein condensation at low tem- 
perature (because true Bose-Einstein condensation can- 
not occur in two dimensions, either) as also discussed 
in the Appendix for the Ji model. Finally, Xio,v) grows 
as ^ l/T without the factor N, and numerically much 
smaller than either X(7r,o) or X(7r,7r)- 

The (0, tt) order, which has the corresponding gap pa- 
rameter equal to zero, plays the role of the non-dominant 
susceptibilities in SBMFT. The Monte Carlo calculation 
in the previous section showed that such non-dominant 
susceptibilities, i.e. X(7r,o) and X(o,o) for 2J2/J1 < 1 and 
X(T,7r) and X(o,o) for 2J2/J1 > 1, are also non- divergent, 
going down to a constant value at zero temperature. Al- 
though X(o,o) behaves in the similar way in both treat- 
ments, X(o.7r) behaves quite differently in SBMFT (di- 
verging like l/T) as it does in classical theories (going 
down to a constant value). Even when a larger spin S 
is used in SBMFT, all three staggered susceptibilities di- 
verge at low temperature and it is not obvious if the 
classical MC results can be fully recovered within the 
SBMFT. Regrettably, little is known in the literature 
about the connection of the SBMFT in the large-S' limit 
to the classical results of the same Hamiltonian. While 
this would be an interesting issue in its own right, we 
leave it as a subject of future consideration. 

The findings of the SBMFT study can be summarized 
as follows. The imiform susceptibility shows a mono- 
tonic temperature dependence that can be roughly un- 
derstood as linearity. The mean-field transition tempera- 
ture Ty^'^^'^ plays a role similar to the crossover temper- 
ature Too in the classical theory, and has a larger value 



when 2 J2/ Ji gets larger (but only above the critical value 
rjc). Secondly, the two staggered susceptibilities X(ir,7r) 
and X(7r.o) show the divergent behavior which is consis- 
tent with that of X{-k,tt) in the Ji model. Both (tt, tt) and 
(tt, 0) (quasi) ordering are present in the SBMFT calcula- 
tion for all ratios of 2 J2/ Ji and only the relative weights 
of the corresponding susceptibilities shift with 2J2/J1. 



IV. SUMMARY AND DISCUSSION 

In light of the recent controversy about how best to de- 
scribe the spin dynamics of the parent FeAs compound 
above the magnetic ordering temperature, we carried out 
a detailed magnetic susceptibility calculation of the J1—J2 
spin Hamiltonian employing classical Monte Carlo and 
Schwinger boson mean field theory techniques over the 
parameter range < 2J2/J1 < 2. In both treatments the 
long-range magnetic ordering at finite temperature is ab- 
sent and a thorough examination of the correlated para- 
magnetic phase becomes possible. In both approaches 
the uniform susceptibility was found to show a largely 
linear temperature dependence that extends up to the 
exchange energy scale T ^ Ji, in agreement with the 
persistence of linear susceptibility up to a few hundred 
Kelvin in FeAs compounds. Our calculation indicates 
that the linear-T uniform susceptibility behavior persists 
over a wider temperature range when 2J2/J1 ratio is fur- 
ther removed from the critical point 2J2/J1 — 1 (classical 
MC) or 2J2/J1 = rjc (SBMFT), and is largely indepen- 
dent of the degree of frustration in the spin model. 

It was recently pointed out that the uniform sus- 
ceptibility feature also exists in a model based on 
the itinerant electron picture^^. To sharply test the 
validity of one model over the other for FeAs, then, one 
would have to go beyond the magnetic susceptibility 
issue to address other physical properties within both 
approaches. It is also desirable to generalize the J1 — J2 
spin Hamiltonian to the doped case where the linear 
magnetic susceptibility was found to persist as well. 
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APPENDIX A: THE SUSCEPTIBILITY 
BEHAVIORS OF THE Ji MODEL AT LOW 
TEMPERATURE 

The antiferromagnetic Heisenberg model with J2 = 
and Ji = 1 is 
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Hs 



Sj, 



(Al) 



out the two points k = (0, 0) and k = (tt, tt) from the 
above sum and treat the rest as an integral, which gives 



has been studied in SBMFT by Auerbach and Arovas-'^'"', 
but it has not been analyzed in detail how the structure 
factors and the susceptibilities behave at low temperature 
for the finite system size N = L x L. In this Appendix, 
we will give a simple explanation why the uniform and 
staggered susceptibilities behave as X(o.()) ~ ^ + BT and 
X(7r,7r) ~ -/V/T at low temperature in SBMFT for finite 
N. ' 

The static structure factor is given by^^ 



-n 



Bk + - 



1 



where the definition of A2k can be found in Eq. (9). 

The energy spectrum is Ek = \/\^ — ^2k' ^^"^ ^^'^ Bose- 
Einstein function is Bk = l/{e^^'' — l). The mini- 



mal value of Ek occurs at A; = (0,0) or fc = (7r,7r), 
and it is convenie nt to int roduce the gap parameter 
6 = min(iJfe) = \/A^"^^4Af. By numerical calculation 
we find at low temperature 



C 



rpl/2 



(A3) 



to an excellent fit and C is about 3.2 {S = 1). 

For q = (0, 0), the uniform structure factor becomes 



(0,0) 



12Bk{Bk + i) = lYl 



1 



smh\pEk/2) 



(A4) 



(0,0) 



1 



N 



l/a 



kdk 



2sinh^(/3(5/2) 

+ T 



■iTTJi/L sinh^(/3A2A;/2) 

l/a 



/I 

2NT NT^ , [L 



(A5) 



where l/a serves as an upper momentum cutoff. Based 
on the relationship of the structure factor and the sus- 
ceptibility in Eq. (18), one has 



(A2) 



X(o,o) 



2 

C2 



T 



In 



L 



(A6) 



Indeed the low-temperature uniform susceptibility is 
given by the form A + BT . 

For q= (tTjTt), the structure factor becomes 



';•+ 

(7r,7r) 



k 



(A7) 



The two isolated fc-points k = (0, 0) and k 
tribute to the structure factor 



4^2 2^2 y2 



(tt, tt) con- 



(A8) 



when we approximate [e^^'^ — 1)^^ ~ T/6. Inserting Eq. 
(A3) yields ^ N"^ and hence the low-temperature 

susceptibility X(t!:,tv) ^ N/T. The contribution from the 



The energy spectrum expanded near k = (0, 0) and k = 

(tt, tt) becomes E^ ~ -^/P + Ap?. It is useful to separate other k points is S'^^ ~ 0{N). 
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